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ABSTRACT 

We explore the possibility that GSC 07396-00759 (spectral type Mle) is a 
widely separated (~2.82', or projected separation ~12,350 AU) companion to 
the "old" (age ~12 Myr) classical T Tauri binary system V4046 Sgr AB, as 
suggested by the proximity and similar space motions of the two systems. If the 
two systems are equidistant and coeval, then GSC 07396-00759, like V4046 Sgr 
AB, must be a spectroscopic binary with nearly equal-mass components, and 
V4046 Sgr must be at least ~8 Myr old. Analysis of a serendipitous Chandra 
X-ray gratings spectrum and light curve as well as XMM-Newton light curves 
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and CCD spectra of GSC 07396-00759 obtained during long exposures targeting 
V4046 Sgr AB reveals a relatively hard (Ty ~ 10 7 K) X-ray spectrum, strong 
flaring, and relatively low-density plasma. These X-ray characteristics of GCS 
07396-00759 are indicative of a high level of coronal activity, consistent with its 
apparent weak-lined T Tauri star status. Interactions between V4046 Sgr AB 
and GCS 07396-00759 when the two systems were more closely bound may be 
responsible for (a) their dissolution ~ 10 6 yr ago, (b) the present tight, circular 
orbit of V4046 Sgr AB, and (c) the persistence of the gaseous circumbinary disk 
still orbiting V4046 Sgr AB. 



1. Introduction 



Thanks to their proximity, individual members of the various nearby (D ~100 pc), young 
(age ~10-70 Myr) ste llar groups identified over the past ~15 yr (for a recent review, see 
Zuckerman et al.ll201ll ) provide readil y accessible examples of, among other things, the late 
evolution of protoplanetary disks (e.g. JAkeson et al.ll2011t iHughes et al.ll201lL an d references 



therein) and the early stages of evolution of exoplane tary systems (IMarois et al.l 120081 ) and 
hierarchical binary star systems (IKastner et al.ll2008af). I n their tests of a method to identify 



such comoving groups of young stars, iTorres et al.l (120061 ) established as a candidate member 
of the f3 Pic Moving Group (/3PMG) the close binary classical T Tauri system V4046 Sgr 
(P ~ 2.4_di " 



Stempe 



method, Torres et al. 



s &: Gahml 120041 . and references therein) . Via the cluster traceback 



( 20081 ) estimated a distance of only 73 pc to V4046 Sgr. Although 
the age of the /5PMG is estimated to be ~12 Myr (IZuckerman fc Sondl2004l . and references 
therein), t he twin members of this fascinating, short-period binary are both evidently still 
accreting (IStempels fc Gahmll2004j) from a relatively large and massive circumbinary disk of 



gas and dust (IKastner et al.ll2008bl ; iRodriguez et al.ll2010i ). 



Torres et al.l (120061 . 120081 ) also identified, as another /3P MG candidate, the star GSC 



07396-00759 (hereafter GSC0739, spectral type Mle or M1.5e; iRiaz et al.ll2006f ). on the basis 
of its proximity to V4046 Sgr (separation 2.82') and the similar radial velocities and simi- 
larly large photospheric Li abundances of the two systems (t he Li A6708 equivalent widths of 



V404 6 Sgr and GSC0739 are 440 and 200 mA, respectively; IRiaz et all 120061 : Ida Silva et al. 



20091 ). Given the low surface density of known /3PMG members, it is reasonable to hypoth- 
esize that GSC0739 is a dis tant (projected sepa ration ~12.4 kAU) companion to the V4046 
Sgr binary, as speculated by lTorres et al.l (120061 ) . In this respect the relationship of GSC0739 
to V4046 Sgr may be very much like that of 2M1235-39 to the well-studied TW Hya As- 
sociation (TWA) system HR 4796AB (IKastner et al.ll2008al ). whose primary is a young (/? 
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Pic-li ke) A star with debri s disk, or that of the brown dwarf candidate TWA 28 to TW Hya 



itself ( ITeixeira et al.ll2008l ) 



If GSC0739 and 2M1235-39 (= HR 4796C, the putative tertiary component of HR 4796) 
are distant companions to the better-studied V4046 Sgr and HR 4796, respectively, then their 
present orbital periods are ~10% of the (~10 Myr) ages of these systems. Hence these objects 
may offer unique insight into the orbital evolution (and especially the dissolution) of young 
hierarchical binaries and, perhaps, the likelihood of, and conditions necessary for, formation 
of planets with circumbinary vs. circumstellar orbits. 

Here, we briefly evaluate the available data concerning the space motions of V4046 Sgr 
and GSC0730, and we analyze optical/IR photometry and serendipitous X-ray observations 
of GSC0739 (including new XMM-Newton X-ray data obtained during the course of a coor- 
dinated observing campaign targeting V4046 Sgr; Argiroffi et al. 2011). We use the results 
to assess the likelihood that GSC0739 is indeed a distant companion to V4046 Sgr and to 
better ascertain the nature of the former. Finally, we consider the implications if the two 
systems are (or were) physically bound. 



2. Data, Analysis, and Results 



2.1. Proper motions and radial velocities 

The UCAC3 catalog JZacharias et alihoiol ) proper motions (PMs) listed for V4046 Sgr 
and GSC0739 are (fi a , fi s ) = (+3.3±1.7, -52.0±1.3) mas yr" 1 and (+1.1±2.8, -40.0±13.6) 
mas yr _1 , respectively, indicating that (given the UCAC3 measurement errors) the PMs of 



the two systems are indistin guishable]. The hel iocentric radial velocities (Vheiio) of the two 
systems are also very similar: iTorres et al.l (120061 ) list Vheiio = —5.7 km s _1 for GSC0739, and 
CO radio line measurements y ield V^iio = —6. 2 ±0.2 km s -1 for V4046 Sgr AB (IKastner et al. 
2008bl : iRodriguez ei~aUl20~1 0)rl 



1 The identical PMs listed for the two systems by ITorres et al.l (|2006l their Table 6) both correspond to 
the Tycho catalog PM of V4046 Sgr AB, i.e., (+2.1 ± 2.1, -54.5 ± 2.3) mas yr _1 ; there is no Tycho catalog 
listing for GSC0739. 

2 Recent results from high- resolution optical spectroscopy indicate that the systemic velocity of V4046 
Sgr AB may vary by ~ 0.5 km s _1 (Donati et al. 2011). 
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2.2. 



Optical/IR SED and comparison with pre-MS tracks 



By matching Kurucz model atmospheres (IKuruczl 119931 ) appropriate for late-type main 
sequence stars (i.e., logg = 4.5) to optical/near-IR photometry available in simbae§ for 
GSC0739 (Tabled]), we deduce a photospheric effective temperature 35 00±100 K, consis- 



tent with the M l-1.5 spectral type previously determined for GSC0739 ( [Torres et al.l 12006 



Riaz et al.l 120061 ). From the Kurucz model normalization, and assuming the distance to 
GSC0739 is 73 pc (i.e., that GSC0739 snd V4046 Sgr are equidistant), we determine a 
bolometric luminosity Lboi = 0.15 L & (with an approximate uncertainty, dominated by the 
uncertainty in distance, of ~ 25%). This estimate is consistent with that obtained from the 
J magnitude and V-K color of GSC0739, L ho \ = 0.17 L Q , following the method described in 
Kenyon fc Hartmannl (119951 ). 



We used the foregoing results for GSC0739 and adopted the luminosity and temperature 
results reported for the individual components of V4046 Sgr AB by Donati et al. (2011) to 
place these systems on theoretical pre-MS tracks (Fig. [T]). The comparison illustrates that 
the two systems can be both equidistant and coeval only if GSC0739, like V4046 Sgr AB, 
consists of two components with nearly equal l u mino sities. Indeed, GSC0739 is flagged as 
a possible spectroscopic binary by iTorres et al.l (120061 ). The positions of the GSC0739 and 
V4046 Sgr AB systems with respect to the overlaid model evolutionary tracks indicate ages 
of between ~10 and ~20 Myr, if GSC0739 is a binary. This is consistent with independent 
(kinemati c and pre-MS evolutionar y track) age estimates of ~12 Myr for the (3PMG (see 



review m 



Zuckerman &: Songi 120041 ) . On the other hand, if GSC0739 is single, then Fig. [T] 
indicates that its age is in the range ~3-10 Myr, i.e., somewhat young for a /3PMG member. 



2.3. X-ray light curves and spectra 



An off-axis Chandra/HETGS gratings spectrum of GSC0739 (total exposure 144.6 ks; 
OBSIDs 5422, 6265; PI: Hercze g) was serendipitousl y obtained in 2005 Aug. during an 
observation targeting V4046 Sgr (jGunther et al.l 120061 ). The GSC0739 X-ray spectral data 
are uncontaminated by V4046 Sgr AB, but the spectral resolution and the exposure of the 
GSC0739 source are somewhat compromised by its (~3') off-axis position. We extracted and 
combined the medium-energy gratings (MEG) and high-energy gratings (HEG) spectra for 
the full ~145 ks exposure duration and generated a light curve covering the first (longer) of 
the two exposure segments; the resulting light curve and spectrum are displayed in Fig. [2j 



3 The SIMBAD database is maintained and operated at CDS, Strasbourg, France. 
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The total (HEG+MEG) count rate, averaged over the entire 144.6 ks exposure duration, 
was 0.014 counts s -1 . 

XMM European Photon Imaging Camera (EPIC) data were serendipitously obtained 
for GSC0739 in 2009 Sept., again during an observation (3 exposures of ~100 ks) targeting 
V4046 Sgr (Argiroffi et al. 2011). As the source fell along a bad column of EPIC's pn detector, 
we only analyzed data from the two MOS detectors (the source also fell near the edge of the 
active area of MOSl, but these data were only slightly compromised). We used the XMM 
Scientific Analysis System (SAS0 version 10.0.0) to extract MOSl and MOS2 light curves 
for the full (~5 day) duration of the 3 x 100 ks exposure (top panel of Fig. [3]) as well as 
MOSl and MOS2 CCD spectra and responses for each individual ~100 ks exposure interval 
(bottom panels of Fig. [3]). Calibrations were performed using the current calibration files 
(CCF) from release note 271, 21-Dec-2010. Combined MOS1+MOS2 count rates are listed 
in Table H 

The Chandra and XMM light curves reveal strong X-ray flaring at GSC0739. An impul- 
sive flare with amplitude ~10 times the quiescent flux level and exponential decay is observed 
starting at about 65 ks into the first ~100 ks segment of the Chandra observation (Fig. HI 
top panel). A flare of similar shape (with amplitude ~4 times quiescent) is observed near 
the beginning of the overall ~5-day XMM exposure, and its decay appears to last the entire 
duration of the first ~100 ks exposure segment; a smaller flare (and several small count rate 
spikes) are also seen in the second and third XMM exposure segments (Fig. [31 top). 

Table [2] and Fig. [3] (bottom panels) summarize the results of fits of variable-abundance 
(see below) absorbed two-component thermal plasma models to spectra extracted from 
the three ~100 ks XMM/EPIC (MOS) exposure intervals. A similar model fit to Chan- 
dra/HETGS data (not shown) yields similar results — i.e., a prominent hard component 
with Tx ~ 20 MK, and a weaker soft component with Tx ~ 5-7 MK — despite the dif- 
ferent instrumentation; we conclude that, overall, the plasma physical conditions did not 
change markedly between 2005 and 2009. Nevertheless, the strong short-term variability of 
the source is apparent in the XMM fit results (Table |2]), which indicate a steady decline in 
the emission measure of the hard component and, possibly, a slight decrease in Tx as the 
observation progressed and the strong flare faded (the source may have suffered from photon 
pileup during the first exposure interval). 

Individual line intensities in the Chandra/HETGS data are not particularly well fit by 
the foregoing simple (two-component) model, as expected given the continuum of plasma 
temperatures that is likely present in the X-ray-emitting region. Hence, we performed global 



4 http: / /xmm.esa.int / sas 
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fitting of a broken pow er law emission meas ure distribution mode l (using isig^ and ATOMDB 
v2.Qj), such as used by lRosner et al.l (ll978| ) or lPeres et al.l ( 1200 ll ) for modeling coronal loops, 
to the gratings data. The model is specified by an emission measure normalization, the peak 
temperature (T peak ), the power law slopes below (a) and above the peak temperature, 
and the elemental abundances. The interstellar absorption was assumed to be negligible (a 
safe assumption given the very modest absorption determined from the fits to XMM data; 
Table [2]). To account for the instrumental off-axis PSF, we used a broadening equivalent 
to a Gaussian turbulent velocity of 600 km s -1 . The resulting best fit model is overlaid on 
the Chandra/HETGS spectum in Fig. [2] (bottom panel); the best-fit model parameters are 



T, 



peak 



11.2 MK (0.97 keV), a = 0.2, (3 = -1.8, and an emission measure of 3.3 x 10 52 cm" 



The best-fit abundance ratios (relative to solar) are Ne/O = 1.8, Mg/O = 0.27, Si/O = 0.6, 
and Fe/O = 0.27 (the two-component plasma model fits to the XMM/MOS data also indicate 
an elevated Ne/O ratio and depress ed Fe/O rati o). These Ne/O and Fe/O ratios are typical of 
coronally active stars (see review in iTestalboioh . Analysis of X-ray gratings spectra of V4046 
Sgr AB likewise indicates enhanced Ne/O a nd low Fe/O, althoug h the Ne/O enhancement is 
more extreme (a factor ~5 larger than solar; iGiinther et al.ll2006l ) in the case of this (actively 
accreting) system. 

Although the resolution of the Chandra/HETGS spectrum of GSC0739 suffers from 
degradation due to the off-axis position of the source, it is apparent that the forbidden 
to intercombination line ratio within the He-like Ne IX line complex (at ~13.5 A) is large 
(Fig. [21 inset in bottom panel), indicative of lower-density (~10 10 cm~ 3 ) plasma. Fig. [2] 
(bottom panel) also illustrates the relatively large Ne x to Ne IX line ratios of GSC0739, 
reflecting the dominance of relatively hot (~10 MK) plasma in its emission measure dis- 
tribution (Table [2]). These density- and temperature-sensitive Ne line ratios are consistent 
with coro nal X-ray emission and, hence, a weak-lined T Tauri classification for this source 
(see, e.g., iKastner et al.ll2004l ; iHuenemoerder et al.ll2007l ). Both Ne line spectral diagnostics 
also stand in stark contrast to those of V4046 Sgr AB; the latter displays Ne (and O) line 
ratios indicative of a signific ant emission contribu tion from cooler (~3 MK) plasma at much 



higher density (~10 cm ; IGiinther et ajj|2006j , Argiroffi et al. 2011), as expected if much 



of the X-ray emiss ion (essentially, all of the soft component) is generated in accretion shocks 
(jSacco et al.ll2010l . and references therein). 



5 http://space. mit.edu/cxc/isis/ 
6 http: / /www. atomdb.org/ 
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3. Discussion and Conclusions 

The available proper motion and radial velocity data (§2.1) support the hypothesis that 
V4046 Sgr AB and GSC0739 are comoving systems. Hence - provided GSC0739 is also a 
close binary (§3.2) - the two systems likely constitute a hierarchical multiple and, in the 
following, we refer to GSC0739 as V4046 Sgr C[D]. 



3.1. The nature of V4046 Sgr C[D] 



Based on its lack of Ha emission (Ha equivalent width of +3.1 A; iRiaz et al.l 12006 ). 
V4046 Sgr C[D] appears to be best classified as a weak-lined T Tauri system (wTTS); i.e., 
unlike V4046 Sgr AB, the system shows no evidence for ongoing accretion. The prominent 
hard (TV ~ 10 7 K) X-ray spectral component, strong X-ray flaring, and relatively low-density 
X-ray-emitting plasma of V4046 Sgr C[D] (§2.3) — all of which are indicative of a high level 
of coronal activity — are fully consistent with such a wTTS classification. 



3.2. V4046 Sgr AB and C[D]: constraints on system age and origin 



If the V4046 Sgr AB and C[D] systems are indeed bound, they are only weakly so, given 
their very large projected separation of ~12,350 AU (assuming D = 73 pc). Assuming the 
systems are still bound, the orbital period of V4046 Sgr C[D] with respect to V4046 Sgr AB 
would be of order ~ 10 6 yr, i.e., a significant fraction of their estimated ages. This system 
would therefore appear to be the most loosely bound hierarchical binary system known in 
the /3PMG. However, two similarly loosely bound binaries have been id entified thus far in the 
TWA: TWA 11AB and TWA 11C, projected separation ~13500 AU fcastner et al.l 2008ah 
and TW Hya (TWA 1) and TWA 28, projected separation ~41000 AU jTeixeira et al.ll2008h . 



The separations of these nearby, young hierarchical multiple systems are comparable to those 
of the most loosely bound binaries in the solar neighborhood ( the separation dis tribution of 
nearby field binaries appears to be truncated at ~20,000 AU; IClose et al.lll990h . The (very 
low- mass) TWA 30A an d 30B star/disk systems, which are separated on the sky by ~3400 
AU (ILooper et al.ll2010l ). are also very weakly bound. 



However, the analysis in §2.2 strongly implies that V4046 Sgr AB and the putative V4046 
Sgr C[D] could only be coeval if the latter, like the former, is a nearl y-equal-component s 
binary — consistent with previous suspicions concerning this system ([Torres et al.l 120061 ) . 
Furthermore, if the two systems indeed formed together, then the position of V4046 Sgr 
C[D] in Fig. [1] appears to place a conservative lower limit of ~8 Myr on the age of the V4046 
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Sgr multiple system, while — on the basis of this comparison of V4046 Sgr C[D] with pre-MS 
tracks alone, and given the present uncertainties in the stellar parameters — the system may 
be as old as ~25 Myr. 

Numerical simulations of hierarchical triple star systems demonstrate that the action 
of Kozai oscillati ons, in combination with tida l friction, s hrinks and circularizes t he or bit of 
the inner binary (IFabrycky fc Tremaind 120071 ) . Although iFabrycky fc Tremaind (120071 ) find 
that the fraction of circularized, close (few-day period) binaries with tertiaries as distant as 
V4046 Sgr C[D] should be negligibly small, they also note that ~40% of hierarchical triple 
systems are "lost" from their numerical sample, due to dissolution. Perhaps V4046 Sgr 
AB and C[D] constitute such a (recently dissolved) hierarchical multiple system that, at a 
slightly earlier epoch, resembled closer, young hierarchical binaries (such as the ~8 Myr-old 
TWA systems HD 98800 and Hen 3-600). If so, the dissolution of V4046 Sgr AB and C[D] 
evidently occurred ~10 6 yr ago and left be hind a close binary with a nearly circular orbit 
(V4046 Sgr AB has an eccentricity e < 0.01; Istempels & Gahml 12004 ). 



Given that V4046 Sgr has apparently well ou tlived the "canonical" few-Myr timescale for 
Jovian planet formation (e.g.. ICurrie et al.ll2009l and references therein), one might further 
speculate that the presence of V4046 Sgr C[D] may explain the longevity of the circumbi- 
nary disk around V4046 Sgr AB — just as the hierarchical natures of the young multiple 
systems HD 98800 and Hen 3-600 may explain the persiste nce and observed properties of 



the c i rcumbinary dust disk in each of these systems (see, e.g.. lPrato et al.ll200ll ; iFurlan et al. 
20071 ; Andrews et al.ll2010l and references therein). Although it remains to work out the 



details in all of these examples, in the specific case of V4046 Sgr it is possible that dynamical 
interactions with component (s) C[D] may have inhibited the formation of gas giant and ice 
giant planets in circumbinary orbits around components AB, thereby preserving the gaseous, 
circumbinary disk from which this close binary system is still accreting. 
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Table 1: GSC0739: PHOTOMETRY 11 



B V R I J 


H 


K 


14.14 12.78 12.13 10.64 9.44 (0.02) 


8.77 (0.04) 


8.54 (0.02) 



a) Photometric data compiled from SIMBAD (sources: Tycho and UCAC3 catalogs, and Torres et al. 2006h . 
Photometric uncertainties available only for 2MASS data. 



Table 2: GSC0739: XMM X-RAY SPECTRAL FIT RESULTS' 2 



(obs. 


C b 


N H 


Ti 


EMf 


T 2 


EM^ 


L x 


log-L x /Lboi 


interval) 




(10 19 cm" 2 ) 


(MK) 


(10 52 cm" 3 ) 


(MK) 


(10 52 cm" 3 ) 


(erg cm" 2 s -1 ) 




(1) 


0.24 


0.0-3.0 


5.2-5.6 


3.1-3.3 


21-23 


3.3-3.5 


9.6 x 10 29 


-2.76 


(2) 


0.12 


2.0-8.0 


5.3-5.8 


2.4-2.6 


17-20 


1.1-1.3 


3.8 x 10 29 


-3.17 


(3) 


0.081 


3.0 d 


5.2-5.6 


2.8-3.1 


14-27 


0.38-0.6 


2.9 x 10 29 


-3.28 



NOTES: a) Fit to absorbed two-component, variable- abundance (APED) plasma models, with parameter 
ranges corresponding to 90% confidence intervals, b) Combined MOS1+MOS2 count rate, c) Plasma 
component emission measure, d) Parameter fixed during fitting. 
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Fig. 1. — The HR diagram positions of V4046 Sgr AB (circles; luminosity and temperature 
data from Dona t i et a l. 2011) and GS C0739 (squares; see §2.2) overlaid on pre-MS tracks 
from ISiess et al.l (120001 left panel) and iBaraffe et al.l (119981 right panel) . The empty square 
indicates the observed position of GSC0739; the filled square indicates its position assuming 
it is a binary system composed of a pair of identical stars. In each panel, the evolutionary 
tracks (solid lines) correspond to masses ranging from 0.2 to 1.0 M & (in intervals of 0.1 M & , 
from right to left), while the isochrones (dashed lines) correspond to ages of 5, 10, 15, 25 
and 40 Myr (from top to bottom). 
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Fig. 2. — Top: Chandra/HETGS light curve of GSC0739, covering the longer (~100 ks) 
continuous observing segment of the two-segment 145 ks exposure. Black: source; grey: 
background. Bottom: Chandra/HETGS spectrum of GSC0739 (black, with grey error bars) 
overlaid with best-fit APED plasma model (red) in which the emission measure distribution 
is assumed to follow a power law (see text). The lower panel shows the residuals of the fit. 
The inset is a blowup of the 12.0-13.9 A region, illustrating the strong emission lines of Ne 
x and Ne IX as well as weak lines of highly ionized Fe. 
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Fig. 3.— Top: XMM/EPIC MOS light curve of GSC0739. Binsize is 1500 s. Bottom: MOS 
spectra corresponding to the three ~100 ks observing intervals illustrated in the light curve 
in the top panel, with best-fit two- component thermal plasma models overlaid. In all panels, 
red points indicate MOS1 data and blue points indicate MOS 2 data. 



